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A SECOND-ORDER SHOCK-EXPANSION METHOD APPLICABLE TO BODIES OF REVOLUTION 

NEAR ZERO LIFT 1 


By Clarence A. Svvertson and David II. Dennis 


SUMMARY 

A second-order shock-exp msion method applicable to bodies 
of revolution near zero lijt »s developed. Expressions defining 
the pressures on noninclined bodies - y derived by the use of 
characteristics theory in combination with properties of the flow 
predicted by the generalized shock-expanlion method. This 
result is extended to inclined bodies to obtain impressions for the 
normal-force and pitching -moment derivatives at zero angle of 
attack . The method is intended for application under iftndi-^ 
tions between the ranges of applicability of the second-order 
potential theory and the generalized shock-expansion method — 
namely , when the ratio of free-stream Mach number to nose 
fineness ratio is in the neighborhood of 1 . 

For noninclined bodies , the pressure distributions predicted 
by the second-order shock-expansion method are compared with 
existing experimental results and with predictions of other 
theories . For inclined bodies , the normal-force derivatives and 
locations of the center of pressure at zero angle of attack pre- 
dicted by the method are compared with experimental results for 
Mach numbers from 3.00 to 6.28 . Fineness ratio 7, 5 , and 3 
cones and tangent ogives were tested alone and with cylindrical 
afterbodies up to 10 diameters long . In general the predic- 
tions of the present method are found to be in good agreement 
with the experimental results . For noninclined bodies , pressure 
distributions predicted with the method are in good agreement 
with existing experimental results and with distributions ob- 
tained with the method of characteristics . For inclined bodies , 
the normal-force derivatives per radian {for normal-force 
coefficients referenced to body base area) are predicted within 
± 0.2 and the locations of the center of pressure are predicted 
within ±0.2 body diameters. On the basis of these results , the 
second-order shock-expansion method appears applicable for 
values of the ratio of free-stream Mach number to nose fineness 
ratio from Ofi to 2. 

INTRODUCTION 

The flow about bodies traveling at high supersonic speeds 
was investigated by Eggers and Savin (ref. 1). They found 
that under specified conditions such flows could be consid- 
ered as locally two-dimensional and could be treated by a 
generalized shock-expansion method. This method has 
been applied to nonlifting and lifting bodies of revolution in 


references 1 through 4. It was found that the generalized 
shock-expansion method accurately predicted the flow about 
pointed bodies of revolution when the hypersonic similarity 
parameter (ratio of Mach number to body fineness ratio) 
was greater than about 1. This method is, therefore, par- 
ticularly useful in the treatment of flows about bodies travel- 
ing at relatively large Mach numbers. At lower speeds, the 
second-order potential theory of Van Dyke (ref. 5) has been 
widely used. (See, also, his hybrid theory for slightly in- 
cf/m'd bodies, ref. 6.) The amylkation of this tl^oiy to 
bodies traveling at urge Mach numbers is often limited, 
however, by the restriction that the maximum slope of the 
body must be somewhat less than the slope of a free-stream 
Mach line. 

The ranges of applicability of the generalized shock- 
expansion method and the second-order potential theory do 
not always overlap, and there remain, therefore, flows at 
certain combinations of Mach number and body shape which 
cannot be treated by either theory. Normally, these in- 
termediate flows are encountered when the hypersonic 
similarity parameter based on nose fineness ratio is in the 
neighborhood of 1 . Since this is a range of practical inter- 
est, additional theoretical methods are needed. Some of 
this need has been fulfilled by the hypersonic small-dis- 
turbance theory (refs. 7 and 8). In its present state of 
development, however, this theory lias application only to 
limited classes of noninclined bodies of revolution. For 
example, due to the series form used to represent the pressure 
distribution, it cannot be applied to the nose-cylinder 
combinations commonly employed for missile bodies. In 
large part, then, the need for a theory applicable at values of 
the hypersonic similarity parameter near 1 still remains. 

The present report develops a theory intended to fulfill 
this need. It is developed by an iteration procedure which 
employs the generalized method of references 1 through 4 as 
the first approximation. This theory is therefore called a 
second-order shock-expansion method. Expressions are 
derived which define the pressures on noninclined bodies of 
revolution. Expressions are also obtained for the normal- 
force and pitching-moment derivatives at zero angle of 
attack. Predictions of the method are compared with those 
of other theories and with experimental results. 


* Supersedes NACA TN 3527 by Clarence A. Syvertson and David H. Dennis, 1956. 
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SYMBOLS 

body cross-sectional area 
body base area 

function defined by equation (B15) 

characteristic coordinates 

. . _ . normal-force 

normal-force coefficient, j 

Qco 

pitching-moment moment about body vertex 
coefficient, " q n Agd 

pressure coefficient, — - " 

body diameter 
entropy 

fineness ratio (Fineness ratio of the nose section 

is/*.) 

body length 
Mach number 
static pressure 
total pressure 
dynamic pressure 

rectangular coordinates (streamline direction and 
normal to streamline direction, respectively) 
cylindrical coordinates (x measured from vertex 
of bo4y and <p from windward meridian) 
en?ter-of-pjressure position (measured from body 
vertex) 

angle of attack, radiates 
function defined by equation (12) 
ratio of specific heats (1.400 for air) 
flow deflection angle 
function defined by equation (9) 
loading (defined by eq. (14)) 
function defined by equation (B16) 

Mach angle (arc sine 1/A/) 

Prandtl-Meyer expansion angle 
shock-wave angle 

function defined by equation (13) 

ratio of cross-sectional area of streamtube to that 

at M= 1 (see eq. (B17)) 

SUBSCRIPTS 

free-stream conditions 

conditions evaluated at various points in flow field 
afterbody 

quantities evaluated for cone tangent to the body 
quantities evaluated by generalized shock-ex- 
pansion method 

quantities evaluated at vertex of body 
quantities evaluated along downstream face of 
shock wave 

quantities evaluated for cones tangent to body 
vertex 

quantities evaluated for cones tangent to body 
at station x 


THEORY 

In the present development of a second-order shock- 
expansion method, attention will be restricted to bodies of 
revolution. It is recognized, however, that the procedure 
used herein may, in the future, find application to other three- 
dimensional shapes. 


The present method is a refinement of the generalized 
shock-expansion method of references 1 through 4. On the 
surface of a body of revolution, immediately behind a corner, 
the generalized method represents a first-order solution for 
the flow and the present method gives the second-order 
solution (see Appendix A). Before proceeding, therefore, it 
is well to orient the present analysis with ;i review of the 
approximations involved in the treatment of the flow about 
bodies of revolution with the generalized method. These 
approximations may be listed as follows (see, e. g., ref. 4): 
(1) Disturbances incident on an oblique shock wave are 
largely absorbed therein, and hence, reflected disturbances 
are negligible; (2) the flow appears locally two-dimensional; 
(3) surface streamlines may be taken as meridian lines. In 
the intermediate range of supersonic speeds of interest he-e, 
the first approximation is particularly well justified (see, ref. 
9), and it will not be considered further. As a consequence 
of the second approximation, a so it km given by the general- 
ized method satisfies the continuity equation only approxi- 
mately. 2 Although th-_. continuity equation does not appear 
explicitly in the following analysis, it is this approximation 
that is refined by the present method. The third approxi- 
mation applies Only to inclined bodies of revolution, since 
for noninclined bodies, meridian lines are, of course, the exact 
streSmJines. In the present investigation, attention will be 
restricted to bodies near zero lift. Under this restriction to 
infinitesimal angles of attack, an analysis has shown that the 
deviation of true streamlines from meridian lines has negli- 
gible. effect on surface pressures. In the following develop- 
ment, therefore, the use of meridian lines as streamlines will 
be retained. 

NONLIFTING BODIES 

The generalized shock-expansion method was developed 
for nonlifting bodies of revolution from the method of char- 
acteristics (ref. 2). This development may be summarized 
with the aid of the equation for the streamwise pressure 
gradient (see, e. g., refs. 2 and 9). 

bp 2yp bd _ 1 bp , , 

bs sin 2/x bs cos n bC\ 

In the generalized method the pressure is considered constant 
along first-family Mach lines (refs. 1 and 4). As a conse- 
quence, the right-hand member of equation (1) is approxi- 
mated by zero, and the equation can be integrated to yield 
the well-known Prandtl-Meyer relation. The objective of 
the present analysis is to refine this approximation to the 
right-hand member of equation (1). To this end, consider 
the flow about a body of revolution which has a pointed nose 
and over which the flow is everywhere supersonic. The 
problem will be simplified by approximating the profile of 
the body with a series of tangents to the original contour 
(see fig. 1). It might be noted that Ferrari (ref. 10) sug- 
gested a similar scheme with a body whose profile was made 
up of chord lines joining points on the original contour. 
While either approximation is permissible, the tangent body 
was selected here so that the conical flow at the vertex will 


2 In the treatment of two-dimensional flows, the first approximation is used but continuity 
is exactly satisfied. 
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^--Tangent body 



Figure 1. Approximation of body shape by a tangent body. 


be correct regardless of the degree of approximation used 
downstream of the vertex. 

The generalized method gives the exact change in surface 
piessure around the corners of the tangent body but predicts 
no change along the straight-line elements. ” The present 
problem reduces to the determination of the pressure varia- 
tion along the conical frustums. 

Since along these surfaces dS/d.$= 0, equation (1) has the form 

dp_ 1 dp 

ds cos n dC\ ® 

The problem now is to find an approximation to the right- 
hand member of this equation. One approximation is a 
power series in Such an approximation would yield, of 
course, the surface pressure as a function of s, and thus, the 
power series could be written equally well in terms of pres- 
sure, p. For reasons which will become apparent later in 
the analysis, a combination of the two series will prove 
advantageous. Thus, the pressure gradient will be written 

^“ = (^ > o+p)(S , o+S i &* + S 2 6* 2 + . . .) 

where P 0 and S n are constants for a particular frustum 
element. This equation is readily integrated to yield an 
exponential variation in pressure 3 

p=_P 0+ c e ( sw +«f+<+--) (3) 

where C is the constant of integration. Conditions must 
now be sought to evaluate the constants in this equation. 
For this purpose, the flow over a single conical-frustum 
element, such as shown in figured, will be considered. The 



Figure 2. — Flow about^a frustum element. 

3 It might be noted that Ehret, Rossow, and Stevens (ref. 11) found that pressures on ogives 
correlated according to the hypersonic similarity law could bo represented approximately by 
an exponential function of distance. 


subscript 1 will be employed to designate conditions imme- 
diately upstream of the corner preceding the frustum; the 
subscript 2, conditions immediately downstream of the 
corner; the subscript 3, conditions at the downstream end of 
the frustum. The first boundary condition to be applied is 
relatively obvious, and it is 

p=p 2 at s=s 2 
Thus equation (3) becomes 


P = — Po+(P2+Po) 


[s 0 (,_ J2)+Si ^Vs,^h 


] 


( 4 ) 


It is also apparent that, if the conical element shown in 
figure 2 were extended indefinitely downstream rather than 
ending at point 3 (see dashed line in fig. 2), then the pressure 
on this extended surface would approach some limiting value. 
If the element were considered to be infinitely long, then the 
only effect the flow upstream of point 2 will have on the flow 
at infinity is to produce an infinitesimally thin lawyer on the 
surface across which the entropy varies. It is readily 
demonstrated, however, that there is no change in pressure 
across this layer and that the flow outside the layer is conical. 
Thus, a second boundary condition is 


V 


Pc 


as .v-- > oo 


where p c is the pressure on a cone tangent to the original 
body at the same point as the frustum element. This 
condition gives the result 


?o=~p c (5) 

provided, of course, the constants, S n , are negative. These 
constants, Sn, are now all that remain to be determined and 
they may be evaluated in terms of the various derivatives of 
pressure evaluated at point 2. Thus, 


$>= 




\d.s‘ / 2 
Pc~P2 


Pc— Vi 


-So 2 


So— 1 


7&p\ 


2L Pc- Pi 


fS 0 3 +3S„S t 


( 6 ) 


at s=s t and so fortli for additional terms. 

At this point, it is well to comment on the choice of the 
series used to approximate the pressure gradient. For 
example, had the series been written entirely in terms of s, 
then a simple Taylor series for the pressure would have 
resulted. While the present approximation can be reduced 
to a Taylor series for small (s— s 2 ), it has added features which 
are desirable. In particular, it is apparent that, within the 
framework of the approximation involved in the use of 
conical frustums rather than a curved body, the boundary 
conditions in pressure at both s=s 2 and s-*co have been met 
exactly. It is primarily because of the particular choice of 
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series which led to the exponential form that it was possible 
to employ the boundary condition at s— > °° . The meeting of 
this condition is relatively important from a practical view- 
point since it serves to limit the possible error involved in 
the use of a series approximation. For example, had a 
Taylor series been used and all boundary conditions applied 
at s=s 2 , then the error for a given approximation could 
increase indefinitely with distance from the corner. With 
the approximation chosen, however, the error must be 
limited since at large s } the solution again becomes exact. 
This added accuracy is not obtained without incurring 
certain restrictions, however; since, as previously noted, the 
coefficients S n must be negative. This condition must be 
satisfied for the present results to be applicable. 

In the following analysis, only the first term, S 0 , in the 
series will be retained, and thus only the first derivative of 
the pressure at point 2 will be required (see eq. (6)). As 
found in Appendix B, an expression defining this gradient is 



»2 



sin 5 2 — sin 5 2 )+ 


> 


B 2 £h 
B\ £ 2-2 



(7) 


where 


B= 


ypM 2 

2(AF— T) 


and Q is the one-dimensional area ratio or 




M 2 


y±l 

2(7-0 


7+ 1 


2 


With the application of these results, the pressure along the 
conical frustum is given by 

p=Pc—(Pc—p2)e~’' ( 8 ) 


where 



x— x 2 

(Vc — P 2 ) C0S 5 2 


(9) 


may be determined from equation (7) since, for the first 
corner, (dp/ds)i is zero. The tangent-cone pressure, p ct may 
be obtained from reference 12 or figure 3. With the various 
factors in equation (8) thus evaluated, the pressure at the 
tangent point and at point 3 (see fig. 2) can be calculated. 
In like manner, the pressure gradient at point 3 can be 
determined by differentiation of equation (8), or 


/dp\ / Pc—P* \ / d P\ 
\d.s/ a \Pc— V%) \ds */2 


( 10 ) 


With the pressure and pressure gradient at point 3 known 
(the Mach number may also be calculated from the pressure 
in the usual manner), the factors in equation (8) may be 
determined for the next element. This process is, of course, 
repeated for each element of the tangent body. 

The procedure just described is not difficult to apply; 
however, further simplification can be obtained by the use 
of a “two-step” tangent body. This body consists of a 
cone tangent to the original body at the vertex and a conical 
surface tangent to the body at the station where the pressure 
is to be calculated. With this two-step body, the second 
surface is a variable depending on the station in question on 
the original body. For this approximation, equation (8) 
becomes 


where 


p=p c -(l>c-p,)e ** 

(ID 

o x sin h v — r cos 8 V 
P r cos 5— x sin 5 

(121 

B a / S2„ sin 5 \ 
~(Pc— P,)\Q, sin bj 

(13) 


The subscript, s, denotes quantities at the station on the 
body as evaluated by the generalized shock-expansion 
method. With equation (11) it is possible to obtain, very 
rapidly, a first approximation to the pressure distribution. 

The second-order shock-expansion method has been de- 
veloped to predict the pressures on a noninclined body of 
revolution. In the following section this method will be 
extended to lifting bodies. 


It is apparent that, in order to apply equation (8), the 
pressure (and Mach number) on the surface of noninclined 
cones must be known. These quantities may be determined 
from the results of reference 2 or reference 12. For con- 
venience, the surves shown in figure 3 have been plotted 
from the results of reference 12. 

By application of equation (8) on successive elements, 
the pressure distribution on the tangent bod} can be de- 
termined. In particular, the pressure at each of the points 
of tangenev may be calculated and applied to the original 
body. The procedure is as follows: First, the elements 
of the tangent body are selected and the coordinates {x,r) of 
each corner determined. The first element is tangent to 
the body at the vertex, and the flow over this element is 
thus conical. For the first corner, then, the pressure, p tJ 
and the Mach number, Mi (see fig. 2), are the same as at 
the vertex of the original body. The pressure, p 2} and the 
Mach number, M 2} may then be determined with the 
Prandtl-Meyer equations. The pressure gradient, (djp/ds) 2 , 


LIFTING BODIES 

For inclined bodies of revolution, a second-order shock- 
expansion method would involve not only a revised expression 
for the pressures, but, in addition, a revised approximation 
to the shape of the surface streamlines. It is recalled from 
the results of Eggers and Savin (ref. 1) that, according to the 
generalized method, surface streamlines may be approxi- 
mated by geodesics. For bodies of revolution, the pertinent 
geodesics are simply meridian lines (refs. 1 and 3). While 
this result is exact for noninclined bodies of revolution, it is 
only an approximation in the case of inclined bodies. A 
refined approximation corresponding to a second-order 
method undoubtedly could be obtained by graphical integra- 
tion of the momentum equations employing the pressure 
distribution given by the generalized method. However, 
it seems at present that this procedure would involve 
extensive calculations. If attention is restricted to bodies 
near a=0, it can be demonstrated that the deviation of the 
true streamlines from the meridian lines will not influence 
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surface pressures. The approximation of meridian lines as 
streamlines can, in effect, be retained and relatively simple 
expressions can be obtained, therefore, for the initial slopes 
of the normal-force and pitching-moment curves. To this 
end, the expression for the normal-force derivative can be 
written 


considered as dependent on angle of attack. By differentia- 
tion of equation (8), there is obtained 

(16) 


dCtf 2ir r l 

doc A B Jo 


A rdx 


(14) 


where the subscript, a=0, has been omitted for simplicity of 
notation and where A is the nondimensional loading on a thin 
disk normal to the body axis and having unit radius. This 
loading A is given by the equation 


, 2 r* d(p/p m ) 

yM^Tj o da 


cos <pd<p 


05) 


The problem then is to evaluate d{pip m )jda. The develop- 
ment given previously which led to equation (3) also applies 
to bodies at infinitesimal angles of attack. Equation (8) 
also applies; however, the variables in this equation must be 


This equation must satisfy the condition = d{p 2 lp m ) 

da da 

at t}= 0 (i. e., x=x 2 ). By the application of this condition to 
equation (16), the last term (involving dy/da) is eliminated. 
The term, d(p 2 /p a> )/da may be evaluated with the aid of the 
Prandtl-Meyer equation 

d(p 2 !pj djpJpA Pi JL_ d{p t ,) -! p 2 d(Pt 2 ) , 

da da p m p h da y p m p h da 1 

Ferri (ref. 13) has shown that the entropy (and hence the 
total pressure, p t ) on the surface of an inclined cone is con- 
stant (independent of <p). When equations (15), (16), and 
(17) are combined, then, the integrals of the terms involving 

dptjdoc and dp t Jda will be zero (since f cos <p d<p= 0). 
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Equation (15) may therefore be written 


A 


-yMj-L _ 


(l- 


d(p c !pj,\i dipxjpj 


da 




X L da 


1 


cos <pd<p (18) 


The only terms in equation (18) that are functions of (p are 
d(pjp m )/da and d(p 1 /p a0 )/da. These two terms may be evalu- 
ated in terms of the normal-force derivative of the tangent 
dC N I 


cone 


1 da 


t and in terms of A L . After performance of the 


necessary manipulations, there is obtained 

A = ( 1 — e~ n ) (tan +^ _n A 1 

d a \tCX *1 


(19) 


It is apparent from equation (19) that dC N /da for cones must 
be known before the loading, A, can be evaluated. Fortu- 
nately, results for cones are available from reference 14 and 
have been plotted for convenience in figure 4. The loading, 
A, may thus be calculated in the same maimer as the zero- 
lift pressures. In this case, Aj for the first corner is simply 

(tan 5 r ) -j v • 
da tcv 



As before, a first approximation to A can be obtained with 
the two-step body. This approximation gives 


A= (!—«-«) (tan 5) 


d_C» 

da 


+ 




In Appendix C, it is shown that equation (20) leads to very 
simple results for certain common body shapes. 

With the loading, A, known, the normal-force derivative 
may be evaluated by integration of equation (14). In like 
manner, the pitching-moment derivative can be determined 
from the equation 


dC m= - 2tt 
da A b d< 


l 

A rx dx 


( 21 ) 


where the contribution to the pitching moment of the varia- 
tion in local axial forces with angle of attack has been neg- 
lected since it is small for slender bodies (see ref. 15). 

A second-order shock-expansion method for bodies of 
revolution has been developed to predict the pressure dis- 
tribution and the normal-force and pitching-moment deri- 
vatives at a~0. The results are relatively simple in form 
and may be applied to a given body with only a moderate 
amount of computations required. Simplified expressions 
based on an additional approximation have also been pre- 
sented which further reduce the amount of work required. 
It should be noted, however, that open-nosed bodies and 
pointed bodies which produce shock waves other than the 
one at the vertex require special forms of the method. The 
necessary equations for these cases are contained in Appendix 
B. In addition, there are several restrictions on the present 
method which should be mentioned. First, it is apparent 
that if the exponential variation of the pressures is to be 
valid, then the pressure gradient just downstream of the 
corner must have the same sign as the pressure difference 
p c -p 2 . This condition, which derives from the previously 
noted restriction that the constants, S nj in equation (4) be 
negative, is given by 0 in the general case and by 0 
for the simplified method. There is an additional restriction 
on the simplified method, and that is that the two-step 
bodies must be real bodies (i. e., the intersection of the two 
tangent lines must not occur at negative values of x or r). 
This condition is given by /3>0. When or /3^=0, 

all equations reduce to those given by the generalized 
shock-expansion method . 

It remains, of course, to determine the accuracy of the 
second-order shock-expansion method and to define its range 
of applicability. There are sufficient data available, both 
from experiment and from characteristic solutions, with 
which the predictions of the method for zero-lift pressure 
distributions can be compared. However for the case of 
lifting bodies, sufficient data are not available, and for this 
reason, the experiments next discussed were conducted. 

EXPERIMENT 

An experimental program was conducted to determine the 
initial slopes of the normal-force curves and the centers of 
pressure for a series of nose-cylinder combinations. The 
tests were designed, of course, to permit a check on the 
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accuracy of the predictions of the second-order shock- 
expansion method just developed. It is recalled that the 
method is intended for application at values of the hyper- 
sonic similarity parameter, MJf n in the neighborhood of 1 . 
The tests cover a range of MJf n from 0.43 to 2.09. 

APPARATUS AND TESTS 

The tests were conducted in the Ames 10- by 14-inch 
supersonic wind tunnel at Mach numbers of 3.00, 4.24, 5.05, 
and 6.28. For a detailed description of this wind tunnel 
and its aerodynamic characteristics, see reference 16. Nor- 
mal forces and pitching moments for the test models were 
measured with a strain-gage balance. The balance consisted 
of a model support sting on which the moments were meas- 
ured at four points. From these four measurements, the 
normal forces and centers of pressure were determined and 
checked. Measurements w r ere made at nine angles of attack 
from —2° to +4° at each test Mach number. At each 
angle of attack, the values of x/d and C N /a were calculated. 
These values were plotted as a function of angle of attack, 
and the intercepts at a— 0 of the resulting curves gave the 
values of dC N /da and x/d at a — 0. 

Wind-tunnel calibration data (see ref. 16) were employed 
in combination with stagnation-pressure measurements to 
obtain the stream dynamic pressures. Reynolds numbers 
based on the maximum diameter of the models were 


Mach number 

Reynolds number 
million 

3.00 

0. 79 

1 4.24 

.72 

6. 05 

.35 

0.28 

. 15 




MODELS 

Cones and circular-arc tangent ogives of fineness ratios 
7, 5, and 3 were tested alone and with cylindrical afterbodies 
having lengths of 2, 4, 6, and 10 diameters. The models 
were made of polished steel and each had a base diameter of 
1 inch. 

ACCURACY OF TEST RESULTS 

Stream Mach numbers in the region of the test bodies did 
not vary more than ±0.03 from the mean values at Mach 
numbers up to 5.05. A maximum variation of ±0.05 
existed at the highest test Mach number of 6.28. 

The accuracy of the test results is influenced by uncertain- 
ties in the measurement of moments and in the determination 
of the stream dynamic pressure and angle of attack. These 
uncertainties resulted in estimated maximum errors in the 
normal-force derivatives and centers of pressure as shown 
in the following table: 


A/co 

dCsjda 

xjd 

3. 00 

±0.15 

±0. JO 

4.24 

±.15 

±.10 

5. 05 

±.20 | 

±.15 

0.28 \ 

±.25 

±.20 


It should be noted that, for the most part, the experimental 
results presented herein are in error by less than these esti- 
mates. 


RESULTS AND DISCUSSION 

NONLIFTING BODIES 

The second-order shock-expansion method has been devel- 
oped primarily to treat flows characterized by values of M„/f n 
near unity. Accordingly, the method has been employed to 
obtain the zero-lift pressure distributions at for 

several different body shapes. In this and all subsequent 
applications of the present method to curved bodies, the 
tangent bodies employed were formed by elements tangent 
to the original bodies at stations x/l n = 0, 0.1, 0.2, . . . , 1.0. 
The tangent-body approximation is required only if the body 
profile is curved since for cone-cvlinders, and for the cylin- 
drical section of any nose-cylinder combination, the present 
method yields results in closed form. 

The results of these calculations are shown in figure 5 along 
with distributions obtained with the generalized shock- 
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(a) Cone-cylinder, M n = 3.00, 

f n = 3 

(b) Sears-Haack body, M m = 3,00, 
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(c) Ogive-cvlinder, M„ = 3.00, 

fn = 3 

Figure 5. — Surface-pressure distributions for various bodies at 
M<oifn~ 1 and a = 0. 
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expansion method (ref. 2). Distributions obtained with the 
method of characteristics (refs. 11, 17, and 18), which are 
considered to be exact, are also shown. It is apparent in 
figure 5 that the present method provides an improvement 
over the generalized method. The differences in the distri- 
butions obtained with the present method and those obtained 
with the method of characteristics are almost indiscernible. 

In figure 5 (c), comparison is also made with the predic- 
tions of the hypersonic small-disturbance theory (ref. 8). 
The curve shown was calculated by three terms of a power 
series representation of the pressure distribution. As noted 
in reference 8, additional terms will be required before this 
method will accurately predict the pressures on an ogive. 
The small-disturbance theory does have a certain advantage 
in simplicity for, if the coefficients of the series expansion 
are known, the pressure distribution can be calculated very 
easily. This advantage is partially offset by the restriction 
that the series method requires the body profile to have con- 
tinuous derivatives up to the same order as the number of 
terms used in the series. With this restriction, the theory 
cannot be applied beyond the nose-cylinder juncture of the 
body (fig. 5 (c)). 

To investigate the accuracy of the present method at 
values of if n other than 1, the comparisons shown in 


figure 6 have been made. Here, the predictions of the pres- 
ent method and those of the generalized method are com- 
pared with experimental results for fineness ratio 3 and 5 
tangent ogives at Mach numbers of 3.00, 4.24, and 5.05. 
The values of jf n range from O.fiO to 1.68. The experi- 
mental results were taken from reference 3. For all cases 
shown, the predictions of the present method are within the 
accuracy of the experimental data. It is also apparent that 
the predictions of the present method tend to approach 
those of the generalized method as M m jj n becomes appreci- 
ably greater than 1. At MJf n — 1.68, for example, the pre- 
dictions of the two methods differ only slightly. 

In figure 6, comparison is also made with t he second-order 
potential theory (ref. 5) for conditions where this theory is 
applicable (i. e , M ca /f n = 0.60 and 0.85). It is somewhat 
surprising that the present method is as accurate as the 
second-order potential theory even at the relatively low 
value of MJf n of 0.60. 

The results presented in figures 5 and 6 indicate that the 
present method fulfills its intended purpose by providing an 
est imate of the pressures on non inclined bodies of revolut ion 
for values of M n jf n near 1. At values of M^jf n as low as 
0.60 the present method provides results comparable in 
accuracy with those obtained with the second-order potential 


o Experiment (ref 3) 
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Figure 6. — Surface-pressure distributions for noninclined tangent ogives. 
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theory. At values of M n /f„ approaching 2, the predictions 
of the present method and those of the generalized shock- 
expansion method differ only slightly. It remains now to 
investigate the applications of the method to inclined bodies. 

LIFTING BODIES 

The experimental results obtained in the present tests are 
given in tables I and II. Predictions of various theories 
are also tabulated. These include the predictions of the 
present method (with various approximations), the general- 
ized shock-expansion method (ref. 3), first-order potential 
theory (refs. 6 and 19), Van Dyke’s hybrid potential theory 
(ref. 6), and Newtonian impact theory (see, e. g., ref. 20). 
Solutions with the second-order potential theory employed 
in the application of the hybrid theory were obtained with 
the aid of reference 21. (Additional results obtained with 
the first-order and hybrid potential theories and with 
Newtonian impact theory may be found in reference 22.) 
With the exception of the two potential theories, all theories 
have been applied throughout the entire range of test 
variables. The potential theories cannot be employed, of 
course, if the free-stream Mach angle is less than the body 
semi vertex angle. 

Normal-force derivative. — The experimentally determined 
normal-force derivatives and the predictions of the various 
theories for the bodies tested are shown in figures 7 (a) 
through 7 (f). Curves for the first-order potential theory 
are not shown in figure 7 since, in all except a few cases, t lie 
predictions of this theory did not differ significantly from 
those of the hybrid potential theory (see tables I and II). In 
general, the present method predicts the normal-force 
derivatives at zero angle of attack essentially within the 
accuracy of the data (within about ±0.2 per radian) through- 
out the entire range of test variables. In addition, the 
present method appears to provide the most consistently 
accurate results of all the theories presented in figure 7. 
The accuracy of the method at low values of M n /f n can be 
explained partially by examination of the predictions of 
the method for the limiting case of very slender bodies. In 
this limit, it can be shown from equations (7) and (9) that 
the term, 77, approaches infinity. From equation (19), then, 
the loading, A, may be written 

A — 2 tail 5—2 ^ (22) 

dC 

since— =2 (see fig. 4). With the substitution of this 
equation in equation (14), there is obtained 



This result is, of course, the well-known prediction of slender- 
body theory, which is known to be accurate for slender bodies 
at low supersonic speeds. Thus, the accuracy of the present 
method at low values of MJf n can be attributed, in part, 
to the fact that it reduces to slender-bodv theory in the 
limit. 


From the results given in figure 7, several observations 
can be made concerning the accuracy of other theories. For 
example, it might be expected that the potential theories 
would he more accurate than the other theories when the 
parameter ^ MJ— 1 tan 5 V is appreciably less than 1. For 
the — 7 cone at = 3 (fig. 7 (a)), however, this parameter 

is only 0.20, and yet, for the longer afterbodies, the hybrid 
potential theory is appreciably more in error than the present 
method. As found in references 2 through 4, the generalized 
shock-expansion method gives accurate results when MJj is 
greater than about 1. Caution should be expressed here, 
however, for the significant parameter is truly M m /f and not 
Mo/jn. The results shown in figure 7 indicate that although 
Mco/Zn may be appreciably greater than 1, for cases where 
the afterbody is sufficiently long to reduce Mm/f below 1, 
the predictions of the generalized method may depart 
appreciably from the experimental results. In general, 
impact theory gives acceptable results only for nose sections 
without afterbodies. 

Center of pressure— The experimentally determined 
centers of pressure and predictions of the various theories 
for the bodies tested are presented in figure 8. The present 
method predicts the location of the centers of pressure essen- 
tially within the accuracy of the data (within about ±0.2 
body diameters) throughout the entire range of test variables 
In addition, the present method again provides the most 
consistently accurate results of all the theories presented. 
In general, all observations made previously regarding the 
reliability with which the various theories predict the normal- 
force derivatives can also be made in the case of the centers 
of pressure. 

Ranges of applicability. — Several parameters are useful 
for defining the ranges of applicability of the various theories. 
The ranges of these parameters covered by the present tests 
are shown in the following table: 


*arameter ! 

Range 

i 

A/ co 

3.00 to 6.28 

f 

3 to 17 

t ft 

3 to 7 

fa 

0 to 10 

Af CD if 

0.18 to 2.09 

ai<±A 

0.43 to 2.09 

1 tan 5. ; 

0.20 to 2.12 


The second-order shock-expansion method was found to 
he applicable throughout the ranges of variables shown in 
the table. Both dC N jda and z/d were predicted within 
±0.2. The present tests did not reveal the limits of appli- 
cability of the method. It was indicated, however, that 
the method may apply to relatively low values of MJi n 
(or \MJ— 1 tan 5,), since, in the limit of very slender bodies, 
the method reduces to the well-known slender-bodv theory. 
The upper limit of the method is dictated by the condition 
specified in the development — namely, 77 > 0 (see eq. (8)). 
Calculations hav e revealed that this condition will be 
violated if \ MJ— 1 tan 5 V is appreciably greater than 2.5. 
It should also be noted that the application of the present 
method to boattailed bodies presents a special problem 















Experiment Second -order shock- Generalized shock- Hybrid potential theory Newtonian impact theory 
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situ-e neither p ( nor is defined in this case (see eqs. 

dot |/ cz 

(8) and (20)). In practice, however, it has been found by 
comparison with results given in reference 23 that the use of 

=2 gives reasonable results for bodies 

tcx 

having moderate amounts of boat tail. 

The present tests also reaffirmed the conclusion given in 
references 1 through 4, that the generalized shock-expansion 
method is applicable when MJj is greater than about 1. 
At values of MJj appreciably greater than 1, no significant 1 
differences between the predictions of the generalized and 
second-order methods were found. The ranges of applica- 
bility of these two methods overlap and thus include most 
Hows about pointed bodies of revolution throughout the 
intermediate- and high-supersonic speed ranges. 

Application of the potential theories is, of course, limited 
by the condition that \MJ— 1 tan must be less than 1. 
Even at the lowest values of 1 tan <5 r covered by the 

present tests, however, neither the first-order nor the hybrid 
potential theory wa? found to provide consistently accurate 
predictions of dCJda or x/d. The calculations performed 
also revealed no significant differen ces in the predictions of 
the two theories at values of yMJ— 1 tan <5, less than 
about 0.7 


APPROXIMATIONS OF THE PRESENT METHOD 

As noted in the development of the present method, a 
simplified solution for bodies with curved profiles can be 
obtained by the use of a two-step tangent body. This 
approximation has been applied to the ogive-cylinders of the 
present tests. By the use of additional approximations to 
the loading, A, the simplified solutions for dC N jda and x/d 
can be obtained in closed form as discussed in Appendix C. 
Examples of the accuracy of the approximate solutions are 
shown in figure 9. While the approximate methods do not 
yield results so consistently accurate as those obtained with 
a more complete solution, the approximate methods may still 
be useful to obtain rapid estimates of dC N jda and x/d. In 
this connection, these quantities can be estimated for ogive- 
cylinders in a few minutes with the aid of the results given in 
Appendix C. 

CONCLUSIONS 

A second-order shock-expansion method applicable to 
bodies of revolution near zero lift has been developed. For 
noninclined bodies, the pressure distributions obtained with 
the method were compared with existing experimental 
results and with the predictions of other theories. For 
inclined bodies, the normal-force derivatives and centers of 
pressure at zero angle of attack determined with the method 
were compared with the predictions of other methods and 
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Figure 9. — Accuracy of approximate solutions for ogive-cylinders (a = 0). 
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with experimental results. Cone- and ogive-cylinders with 
fineness ratios from 3 to 17 were tested at Mach numbers 
from 3.00 to 6.28, corresponding to a range of values of the 
hypersonic similarity parameter based on nose fineness ratio 
(i. e., the ratio of free-stream Mach number to nose fineness 
ratio) from 0.43 to 2.09. These comparisons led to the fol- 
lowing conclusions: 

1. For noninclined bodies, the present method predicts the 
pressure distributions within the accuracy of experimental 
results. At values of the hypersonic similarity parameter 
based on nose fineness ratio as low as 0.6, the present method 
is as accurate as the second-order potential theory. At 
values of the parameter approaching 2, the predictions of 


the present method differ only slightly from those of the 
generalized shock-expansion method. 

2. For inclined bodies, the normal-force derivatives and 
the locations of the center of pressure at zero angle of attack 
predicted with the present method are in good agreement 
with the experimental results throughout the entire range of 
test variables. Within this range, the present method 
yields results more consistently accurate than those of 
other available theories. 

Ames Aekonautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Oct. 12, 1955 


APPENDIX A 

POWER SERIES REPRESENTATION OF FLOW ABOUT BODY OF REVOLUTION 


The accuracy of the present method has been demon- 
strated by comparisons made over h wide range of flow T 
parameters. It is also informative, however, to examine 
briefly the mathematical accuracy of the method. For this 
purpose, the model showm in figure 10 is useful. From the 
vertex to point 1, the body is conical. Between points 1 and 
2, the surface is deflected by a small angle, e. At any point 
downstream of point 1, the physical deviation of the body 
from a conical surface may be given in terms of the angle, e, 
and the distance, As, measured from point 1. Similarly, 
flow parameters at any point downstream of point 1 may be 
expressed in terms of e and As. Before developing such an 
expression, it should be noted that for this model (and within 
the restriction that the flow is everywhere supersonic), the 
present method provides an exact solution for the surface 
flow in several. limits. For example, the present method is 
exact for all values of As when e=0. For As— 0 and As— > o° , 
the method is exact for all values of e. For arbitrary values 
of As and e, of course, the present method is not exact. 
However, the general accuracy of the method (‘an be demon- 
strated by expressing Row parameters in the form of a 
Taylor series in the two independent variables, e and As. 
The dependent variable used to define the flow may be any 
one of several parameters. Pressure and velocity are among 
those most commonly used. In the present analysis, how T - 
ever, the Prandtl-Meyer angle, v , is considered the dependent 
variable. It should be recognized that the value of the 
Prandtl-Meyer angle at a point will define the Mach number, 


A- 

-< 



pressure, velocity, and other such parameters. We have 
then the series 



(As) + 



e+ 




<4!)<,1+ (10, ( *' ) ] + 


(At) 


Each of the derivatives is evaluated at As — e—0. When 
e=0, it is apparent that the flow parameters are constant 
along the surface and independent of s. Therefore, all 
derivatives with respect to s alone are zero. When As=0, 
it is also apparent that (c>^/c><5) x = — 1 and that all higher 
derivatives with respect to 8 alone are zero. We have then 
the problem of evaluating the cross derivatives. The 
second-order cross derivative, (d 2 v/ds d<5), may he evaluated 
with the aid of equation (B14), from which ( dv/ds) 2 may be 
determined; namely, 


/ dA — l/dp\ _ ~ 1 / 12 1 

\ / o A o / 2 2/ 1 ] — 1 \ 1^2 


sin — sin 



(A2) 


It is also apparent that 


/ &P \ ^ d/dA 
\dsb5/i d<5\ds/ 2 


(A3) 


in the limit as Hence, by virtue of equation (A2), 

/ d 2 v \ —cosdi r/ , % , 

( ra ),"^ VW = T [(v ' A, ' ! - lla,,8l) - 11 <A4) 


As noted in Appendix B, equation (B14) is not an exact 
solution for the pressure gradient since a term, 

X 2 f 5 J- / ds 

2 cos ju 2 J 4 A cos fx \bCi/ b 

w r as neglected in the derivation (see eq. (Bll)). In the 
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present analysis, both dpjdCi and the interval of integration, 
(n a — ,s * 4 )/6 (see fig. 11) , are of order e ; lienee, the neglected term 
is of order e 2 . Equation (A2) is accurate to the first order 
in €, therefore, and equation (A4) is exact. With the sub- 
stitution of equation (A4) in equation (Al), and with the 
application of the other results previously noted, there is 
obtained 

*=■= * - 1 tan 5 >) - 1] (AS) (t) + 

0[(A*) 2 (e) , (A«)(c) 2 ] (Ao) 

Tlie generalized shock-expansion method of references 1 
through 4 gives the result that v=v 1 — e. The generalized 
method gives the Prandtl-Meyer angle mathematically 
accurate to the first order of the independent variables e and 
As and, therefore, immediately downstream of the corner, 
gives a first-order solution for the surface flow. The present 
method adds the coefficient of the term involving (As) (e) in 
equation (A5) and, hence, gives the Prandtl-Meyer angle 


mathematically accurate to the second order of the independ- 
ent variables e and As. In general, therefore, immediately 
downstream of the corner the present method gives a second- 
order solution for the surface flow, and therefore, it has been 
termed the second-order shock-expansion method. 

The foregoing analysis considered only expanding flows 
about the corner. If e is positive, then the shock wave 
emanating from the corner must be considered. The result 
obtained is essentially the same, however. For positive e a 
term of 0(e 3 ) must be added to equation (A5) to account for 
the difference between the Rankine-Hugoniot equations and 
the Prandtl-Meyer equations. Alternately, the term, — c, 
in equation (A5) can be replaced with the change in Prandtl- 
Meyer angle between points 1 and 2 as given by the Rankine- 
Hugoniot equations. The second-order term in either case 
is identical, however, as equation (A4) may also be obtained 
by differentiating equation (B21). (It may also be obtained 
by differentiation of the exact pressure-gradient equation, 
eq. (B18).) 


APPENDIX B 

EVALUATION OF PRESSURE GRADIENT DOWNSTREAM OF CORNER OF BODY OF REVOLUTION 


CONVEX CORNER 

Along a streamline in axially symmetric flow the following 
relation holds (see eq. (1)) 


dp ^d<5_ 1 dp __ — X / d<5 _|_sin p sin <5 

ds ds cos jj. d( \ cos p\d(\' r 

From this equation, we may also write 

dp /dp db\ 

d/r ( ‘ OSM fe“ x w 

and 

sin fi sin <5 cos p /dp d8\ 

ds) 


) (Bl) 


dfi 

dCY 


(B2) 


(B:j) 


Consider now the flow in the region of a convex corner on a 
body of revolution as shown in figure 11. Between points 
4 and 5, we may write, 

,r-Mach lines 



from equation (Bl) 

If points 4 and 5 are near to the surface, equation (B4) may 
be approximated by 


dp p 5 ~p 2 Pa— 


f‘ ‘ ($)* 

I Ja x COS p \dCi/ 


IBS) 

Since the flow between points 1 and 2 is strictly of the 
Prandtl-Meyer type, 


r>d V 

J — j/ 2 =5 2 — di 

We may also write from equations (B2) and (B3) 


(BQ) 


P*Z.P2 = 

-J, 

( bp ' 

\ A— C0s 

x 2 

X 2 ’ 

KbC, y 

V x 2 ~ 

1 

s. 

1 

-li 

( bp'' 

v COS Ml 

X, 

X, 1 

KbC t/ 

\ a ~ xr 


x js \ ( \ k * sin M 2 sin 6 2 , cos Mg f/dA . /d5\ T 

!s 4 U >/' * — ; — b ~— LWr’-UAj 6 

(B9) 

. . . / \ „ sin mi sin Si cos mi T/ bp\ ^ / b5\ "1 

s ^ s ‘+{sc,)r > ‘-- -r— -■ TrL(s)r x '(w,> 

(BIO) 

When equations (B6) through (BlO) are substituted into 
equation (B5), there is obtained 
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2 cos 
X 2 


— »=rx-c^(^)* < bu > 


sin fj .2 sin 5g £__sin Mi sin 


If it is assumed that a first approximation to the flow is 
given by the generalized shock-expansion method, then the 
right-hand member of equation (Bll) may be neglected. 
Equation (Bll) may thus be written 


(*P\ /aa\ X 2 r/a\ 
\i>s /2 2 \i>s/ 2 2rcos/z 2 L\^/ 


i 2r cos LVlJ Sl " Si " 5, ““ n * 8 '" { *]+ 

£2 (mi©; -'GO,] < B ' 2 > 

In the limit as the streamline between points 4 and 5 ap- 
proaches the surface, the ratio, a/b, may be evaluated in 
terms of the one-dimensional area ratio 


where 


and, of course 


B-- 


7 pM 2 


U = 


M 


2(M 2 — 1) 
sin 2/x 

1 + ( 7 7 1;) M 2 

~W)~ 


(Y+l) 


2 ( 7 - 1 ) 


(Hi 5) 

(Bl6) 

(Bl 7) 


a / si n ^ 2 \ ft, 
b \sin iij S2 2 


(B13) 


With the substitution of equation (Bl3) into (Bl2), there is 
obtained after combination of terms 


Gf)+= 


/cM\ = B 2 /lli 
\ds/ 2 r \il 2 


sin 5j — sin 5 2 )+ 


2 )‘ 


: I [(§!)+■(£),] <*»> 


Bo 

B 


Equation (B14) represents only an approximate evalua- 
tion of the pressure gradient. More exact evaluations may 
be found in references 10 and 24. These more exact results, 
of course, require numerical or graphical integration. 

CONCAVE CORNER 

In most cases, the tangent bodies used in the application 
of the present method will have convex corners. There is 
a possibility that concave corners may be encountered. 
In the event that the original body does not have sharp 
concave corners, equation (B14) will still suffice since the 
flow along the surface is still isen tropic. However, if the 
original body does have sharp concave corners, then the 
pressure gradient for this case will also be required. This 
result can be obtained in the same way as equation (B 1 4) ; 
however, the shock wave emanating from the^ corner must 
be considered after the manner described in reference 9. 
The expression defining the pressure gradient in this case is 


where 


“? (M). |sin sin s -“ sin <'-« sin SJ+ (aO.{i (M).[ 8J BuT?]+(|r f ) 008 <'-*■>}- 

• Q{(^)R ^] + 1 GO. - (£). - <-« 


(^l) P«M U 2 sin (<t— 5„) cos (< t— 8 u ) 


sin (cr— 5„) cos (<r— 5<j) . , 

( 4 > 

COS 2 (<J—8a) 1 

L sin (cr — 5 W ) cos (<r— 8 U ) 

VY+1/ 

M u 2 sin 2 (<j— 8 U ) J 


(B18) 


(B19) 


and 


Mu ‘ ! ) sin 

f [(7+1) tan (Sg— 8 U ) cos (tr— 3 U )— sin (o— 5 U )] M„ 2 sin 2 (o— 5„)+sin (o— 5 U ) 

\ 1 + [1— 2 sin 2 (cr — 5^) +2 tan (8 d ~ 8 U ) sin (<r—8 u ) cos (cr — 5„)] M u 2 sin 2 (<r— 8 U ) j 


(B20) 


In these equations, <r is the shock- wave angle with respect to 
the body axis, and (dp t /dn) u is the variation of the total 
pressure normal to the surface just upstream of the shock 
wave. The subscript, u y refers to conditions upstream of 
the shock v^ave, and the subscript, d , refers to conditions 
downstream. Equation (B18) represents the exact solution 


in the usual sense. All effects of the interaction between 
shock waves and Mach waves are therefore included. In 
order to be consistent with other parts of this analysis, 
these effects should be neglected. In addition, since equa- 
tion (B18) is intended for application to a tangent body, the 
body curvatures, (d8/ds) u and (d5/d s) d , will be zero. It 
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may also be noted that the first step of the tangent body is 
a cone tangent to the vertex of the original body. For this 
approximation then, there will be a small layer near the 
surface of the tangent body for which (dp t /dn) = 0. With 
these approximations applied to equation (B18), a simplified 
result can be obtained which will suffice for the present 
purposes. 


fdp\ f tan na H 2B d rs\ 
\ds ) d L ' tan(o— 6 d )J r L 


sin(<r— 5 u )sin 5 U _ 


sin(<r — 5 d ) 


sin 5 d + 


]- 


P sin (a i / Pd 3 u )tan /i d (g2i) 

) sin(<7 — 5 d ) J 


For a body with a concave corner, a special form must 
also be used for the loading. Just downstream of the corner 
and before the first convex corner 


A = 


= (l-e-') tan 



(B22) 


and, thereafter, 


A — (1 — e~ v ) tan 5 


dCff 

da 


+^e"*A,+ 

\tcx Ai 



^2 Pi \ 
Xl PcJ 


GA u 


(B23) 


where 


G= 


F 

(Pd/pJ 


{[(7- 


[M u 2 sin 2 {a 6 U ) 1 j 2 


l)M„ 2 sin 2 (<7— 5 U ) + 2]M u 2 sin 2 (<r 


-«»)} 


(B24) 


The corresponding equations for open-nosed bodies of 
revolution are similar. The pressure gradient at the leading 
edge may be determined from equation (B21) with M u =M m , 
5 m = 0, and (dp/d$) u = 0. The loading on the exterior surface 
is given by 


A= (1 — e~ v ) tan 5 


dC N 

da 


+r 2 - e ’Ar 
tcx *1 


-( 


Pi. 

Pco 


X 2 Pl\ 

Xi pj 


JA V 


(B25) 


where A 0 is the loading at the leading edge, or 


A 


4 sin o p cos o v 


(7+1) 1 

sin (<r, — S v ) cos (<r„— 5 P ) 

-( 4 1 

cos 2 (<r„ — SJ - ] 

sin o p cos o p 

\y+V 

' MJ sin 2 <t c J 


(B26) 


and J is defined by 
J= 


(M a 


l) 2 


(pJpjMJ sin 3 o r [(y— l)MJ sin 2 o v + 2] ® ^ 

For bodies with concave corners, and for open-nosed 


bodies, the total pressure is not constant on the surface 
when the bodies are inclined. This variation in surface 
total pressure leads to the term involving A u in equation 
(B23) and the term involving A v in equation (B25). 


APPENDIX C 

EXTENSIONS OF THE APPROXIMATE METHOD 


This analysis is based on the approximate or two-step 
method previously mentioned. The basic equations of this 
method are equations (11), (12), (13), and (20). Before 
proceeding with this analysis, it is convenient to write down 
the expressions for the function (see eq. (12)) for several 
types of bodies. These expressions are presented in the 
following table: 

EXPRESSIONS FOR fi 


Body 

For nose section 

For cylindrical afterbody 

Any body 

x sin <5 ,— r cos 3, 
r cos 3— x sin 5 

2/n sin 5, —cos 3, 

Cone-cylinder 

(Not required) 

* 0-0 
Vi +'*/.* 

Tangent-ogive- 

cylinder 

1 

,+ w ( s K - 0 

1+4/*’ 

Tangent-paraboloid 

cylinder 

>+( ‘- 0 * 


Y /n J +l 

V /»*+ 1 


In general, the equations for the normal-force and pitching- 
moment derivatives may be integrated in two parts — one 
part for the nose section and one for the afterbody. Thus, 
with the loading defined by equation (20) 


dCx 

a. 

1? 

da 

da \ 




(Cl) 


dC m _ 

d(J m 

da 

da 


-S e-OfA [( 1 + G 2 f n ) - (1 + G 2 f n + G 2 f a )e-<hl*} 

nose 


(C2) 


where 


g* a cos 3„ 


(C3) 


and 


6r 2 = 2^ a sin 5„ ^ (C4) 


The additional subscript, a, refers to functions evaluated 
for the afterbody (i. e., 5=0). Thus, from equation (13), 


y(P. a /P«,) M Sa 

4 ' a ~2{i-p,jpj{M !a 2 -i) n Ja 


(C5) 


The terms (?, and G 2 are functions of M m and 5„ alone. 
These functions have been evaluated and the results are 
shown in figure 12. 
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Semivertex angle.S^ ,deg 


(a) Function, G x 

(b) Function, (? z 

Figure 12. — -Curves defining functions G x and Gi (Appendix C). 


For the special case of cone-cylinders, equations (Cl) and 
(C2) represent a closed solution of the general method as 
well. 

By the use of an additional approximation to A, results 
in closed form can also be obtained for ogival nose sections. 
Such an approximation is 


A=(r 3 tan 5 



tcv 


-\-G 4 tan 5 



where 


and 


n _Ps a sin 2(i v 

^7 3 " ' rj ^ 

Pv sin 2fx Sa 


(? 4 — 2(1 — e~* a ) 


(C7) 


(C8) 


When equation (C6) is substituted in equations (14) and (21), 
equations are obtained in closed form for dC N /da and 
dCJda. These equations involve constants which are com- 
plicated functions of the nose angle 8 V (or nose fineness ratio 
f n ). These functions can be expanded in a series in terms of 
8 V ; the leading terms of these series are constants independent 
of 8 V . In view of the approximate nature of this analysis, 
the use of the leading terms will suffice, Thus there is ob- 
tained 


dC N 

da 


ogive 


.+4e.) 


1 dC m 

f n da 


= — — G 4 - 

(3+22(73') 

(dC N \ 

15 4 

15 1 

\ da ) 


(C9) 

(CIO) 


To the accuracy of this analysis, these equations also repre- 
sent the solutions for a tangent paraboloid. These equa- 
tions have been evaluated for a range of Mach numbers 
and nose fineness ratios. The results are presented in figure 
13. It is apparent that with the aid of equations (Cl), (C2), 
(C9), (CIO), and figures 12 and 13, dC N /da and dOJda for 
ogive-cylinders can be evaluated approximately in a few 
minutes. 


da 


2.3 

2.2 

2.1 

2.0 

1.9 

1.8 

1.7 

1.6 




Figure 13. — Curves defining normal-force and pitching-moment 
derivatives for tangent ogives as predicted by approximate closed 
solution of Appendix C. 
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6.28 

0 





2. 23 

2. 11 

1.99 

2. 37 

1.95 

2. 12 


-- 

2 

— 



2.70 

2.47 

2. 36 

3.16 



2. 57 


-- 

4 





2. 97 

2.71 

2. 59 

3.94 


2.78 


-- 

6 

— 



3.13 

2. 86 

2. 74 

4.73 



3.01 


-- 

8 

— 



3. 23 

2. 95 

2.83 

5.52 






10 

— - 

.... 

3.28 

3.01 

2. 89 

6.30 

.... 

3. 31 

Ogive, 

3.00 

0 

.... 


2.19 

2. 06 

1.99 

2.50 

1.86 

2. 28 

A— 3 

-- 

2 

— 



2. 87 

2.60 

2.53 

4.06 



2.98 



4 

— 



3.13 

2.82 

2. 75 

5.63 


3.25 



6 

— 



3.23 

2. 92 

2. 84 

7. 19 



3.20 



8 





3. 27 

2.96 

2.88 

8. 76 




- 

10 

.... 

.... 

3.28 

2.97 

2.90 

10.32 

— 

3. 32 


4.24 

0 





2. 10 

1.98 

1.83 

2. 16 

1.86 

2.09 


-- 

2 

— 



2. 76 

2. 46 

2.31 

3.28 



2. 81 


-- 

4 





3.08 

2. 70 

2.56 

4.40 



3.11 


-- 

6 

— 



3. 24 

2.83 

2.68 

5.53 



3.29 


-- 

8 





3.32 

2.89 

2. 74 

6.65 





- 

10 

.... 

.... 

3. 36 

2.92 

2. 78 

7. 77 

— 

3. 34 


5. 05 

0 

.... 



1.98 

1.88 

1.74 

1.97 

1.86 

2. 01 


-- 

2 





2. 59 

2. 32 

2.17 

2.90 



2. 66 


-- 

4 





2.90 

2. 55 

2.41 

3.82 



2.96 


-- 

6 

— 



3. 05 

2.68 

2.53 

4. 74 


3.16 


-- 

8 

— 



3. 13 

2.75 

2.60 

5. 67 





- 

10 

— - 

.... 

3. 17 

2. 79 

2.64 

6.59 

— 

3.38 


6.28 

0 





1.79 

1.72 

1.65 

1.77 

1. 86 

1.83 


-- 

2 



— 

2. 32 

2.06 

1.99 

2.50 



2.44 


.. 

4 





2.58 

2.26 

2.18 

3.22 



2.61 



6 

— 



2.70 

2.36 

2.29 

3.94 



2.67 


— 

8 

— 



2. 77 

2.42 

2.35 

4. 67 





1 

-- 1 

10 

— 

.... 

2.80 

2.45 

2.38 

5.39 

— 

2.90 
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TABLE II.— CENTERS OF PRESSURE AT ZERO ANGLE OF ATTACK 


Nose 

shape 

A/co 

/. 

First- 

order 

potential 

theory 

Van S 

Dyke 
hybrid 
potential ex 
theory n 

Cone, 

3.00 

0 

4. 67 

4. 67 

/■ = 7 


2 

5. 32 

5. 34 



4 

5.34 

5. 32 



6 

5.00 

5.00 



8 

4.63 

4.63 


-- 

10 

4. 42 

4.41 


4. 24 

0 

4. 67 

4. 67 



2 

5.41 

5.41 



4 

5. 72 

5. 72 



6 

5.69 

5.68 



8 

5.38 

5.38 



10 

5.00 

5.00 


5. 05 

0 

4.67 

4.67 



2 

5.46 

5. 44 



4 

5. 92 

5. 89 



6 

6.04 

6.00 



8 

5.91 

5.87 


- 

10 

5.62 

5.58 


6.28 

0 

4.67 

4. 67 



2 

5. 46 

5. 42 


__ 

4 

6.04 

5.98 



6 

6.38 

6.30 



8 

6.48 

6. 40 


-- 

10 

6.40 

6. 32 

Cone, 

3.00 

0 

3.33 

3.33 

5 

.. 

2 

4.02 

3.98 


.. 

4 

4.09 

4.05 



6 

3.78 

3.75 



8 

3.41 

3. 37 


-- 

10 

3.16 

3.13 


4.24 

0 

3. 33 

3.33 



2 

4.07 

4.05 



4 

4.46 

4.41 



6 

4.45 

4.41 


__ 

8 

4.20 

4. 16 


-- 

10 

3.84 

3.81 


5.05 

0 

3.33 

3.33 



2 

4. 10 

4.04 



4 

4. 59 

4.52 



6 

4.80 

4.72 



8 

4. 74 

4.65 


-- 

10 

4.48 

4. 42 


6.28 

0 

3.33 

3.33 


.. 

2 

4.13 

4.04 


.. 

4 

4. 72 

4.60 



6 

5.14 

5.00 



8 

5. 31 

5. 16 


- 

10 

5.29 

5. 14 

Cone, 

3.00 

0 

2.00 

2.00 

A= 3 

.. 

2 

2. 67 

2.63 



4 

2.85 

2.80 



6 

2. 61 

2.57 



8 

2.25 

2. 20 


-- 

10 

2.00 

1.96 


4.24 

0 

2.00 

2.00 


__ 

2 

2. 74 

2.65 



4 

3.20 

3.08 



6 

3. 34 

3.26 



8 

3.23 

3. 12 


-- 

10 

2.94 

2.84 


5.05 

0 

2.00 

2.00 



2 

2. 77 

2.65 



4 

3.35 

3. 19 


__ 

6 

3. 67 

3.50 


.. 

8 

3.74 

3.56 



10 

3.59 

3. 43 


6.28 

o 




2 





4 





6 




" 

8 

10 
















order 


4. 67 

5.39 

5. 78 
5.99 
6.11 
6. 17 

4. 67 

5.40 

5.88 
6.20 

6.40 
6.53 

4.67 

5.39 

5.89 

6. 25 
6.51 
6.69 

4.67 
5.35 
5.88 
6.29 
6. 59 
6. 81 


General- 
ized 
shock - 
ex pa ns ion 
method 


3.33 
4.02 
4. 43 

4. 66 
4.79 
4.86 

3.33 

4.01 
4.50 
4.85 
5.08 

5. 24 

3.33 

4.00 
4.52 
4.92 
5. 21 
5. 43 

3.33 
3. 95 

4. 48 

4.91 
5.25 

5. 52 

2.00 
2.62 

3.02 
3.27 
3.40 
3.48 

2.00 
2.58 
3.05 
3.38 
3. 61 
3. 77 

2.00 

2. 55 
3. 02 

3. 42 
3. 71 

3.92 

2.00 

2. 49 

2.98 
3.40 

3. 75 

3.98 


4.67 

5. 69 

6.70 
7. 70 

8.71 
9. 72 

4.67 
5. 61 
6.59 
7. 57 
8.55 
9. 55 

4.67 

5. 57 

6. 40 

7. 53 

8. 45 
9 48 

4. 67 
5.50 
6.41 
7. 35 
8.30 

9. 27 


New- 

tonian 

impact 

theory 


Experi- 

ment 


3.33 

4. 29 
5.27 
6. 26 

7. 25 

8. 24 

3.33 

4.21 

5. 14 
6. 10 
7. 07 
8.06 

3.33 

4.15 

5.06 
6.00 

6. 95 

7.92 

3.33 

4.07 

4.93 
5.83 
6. 77 
7. 73 

2.00 

2.85 

3.79 

4. 76 

5. 74 

6. 72 

2.00 
2. 75 
3.63 
4.56 

5.52 
6. 49 

2.00 

2.68 

3.52 
4.43 

5. 37 

6. 33 

2.00 

2. 59 

3. 37 
4.23 
5. 13 
6. 07 


4.67 


4.67 


4.67 


3. 33 


4.75 
5. .56 
5.63 

5.85 

6. 32 

4.01 

5. 55 
6.00 

6. 40 

6.39 

4.71 

5.45 

5.86 
6.28 

6.55 

4. 74 

5. 40 
5.80 
6.31 

6. 74 

3.35 

4.04 

4.46 
4. 53 

4.88 

3. 26 

4. 02 
4. 47 
4.98 

5.39 

3.38 

3. 99 

4. 54 
5.08 

5.25 

3.25 

4. 12 
4.48 
4.95 

5. 32 


2.00 


2.00 


3. 51 

2. 19 
2.68 
3. 18 

3. 35 

3.82 

2. 10 
2. 52 
3.23 
3. 37 

3.82 

2.03 
2. 52 
3. 16 
3.50 

3. 77 


Nose 

shape 

A/co 

/. 

First- 

order 

poten- 

tial 

theory 

Van 

Dyke 

hybrid 

Second-order shock- 
expansion method 

Gener- 
alized 
shock - 
expan- 
sion 
nethod 

New- 

tonian 

mpact 

theory 

Experi- 

ment 

poten- 

tial 

theory 

s 

Ten- 

step 

olution 

Two- 

step 

approx- d 
motion 

Ap- 1 
pen- i 
ix C 

Ogive, 

3.00 

0 

3. 69 

3.66 

3. 73 

3.59 

3.58 

4. 17 

3.17 

3. 79 

h = 7 


2 

3.85 

3. 82 

4. 23 

3.93 

3. 91 

5. 16 


4. 12 



4 

3.62 

3.59 

4. 50 

4. 12 

4. 10 

6. 15 


4. 32 



6 

3.29 

3. 26 

4.63 

4. 23 

4. 21 

7.14 


4. 40 



8 

3.04 

3.00 

4.70 

4. 28 

4.26 

8. 14 




__ 

1ft 

2. 97 

2. 93 

4. 74 

4.31 

4.29 

9. 14 


4. 55 


4.24 

0 

3.80 

3. 76 

3.75 

3.66 

3.63 

4. 11 

3.17 

3.70 



2 

4.26 

4.21 

4.30 

4. 10 

4.06 

5.05 


4.08 



4 

4.34 

4.29 

4.65 

4.40 

4.35 

6. 02 


4.76 



6 

4.20 

4. 15 

4.87 

4.60 

4.55 

6.99 


4.98 



8 

3. 86 

3.80 

5.00 

4.73 

4.68 

7. 97 




__ 

10 

3. 53 

3. 48 

5. 08 

4.81 

4.76 

8.96 


5. 16 


5. 05 

0 

3.85 

3. 77 

3.75 

3. 66 

3. 64 

4.04 

3.17 

3.75 



2 

4.40 

4.31 

4. 34 

4. 15 

4. 13 

4.96 


4.32 



4 

4.61 

4.53 

4.73 

4. 52 

4.50 

5. 91 


4.86 



6 

4. 57 

4. 48 

5. 01 

4. 79 

4.76 

6. 87 

— 

5. 10 



8 

4 34 

4.26 

5.19 

4.98 

4.96 

7. 85 





10 

4.04 

3.95 

5.31 

5. 12 

5.09 

8.83 


5.23 


6.28 

0 

3.88 

3. 74 

3. 72 

3.64 

3. 62 

3.97 

3. 17 

3.65 



2 

4. 50 

4. 37 

4. 31 

4. 14 

4. 14 

4.85 


4.20 



4 

4.86 

4.74 

4. 74 

4.55 

4. 55 

5. 77 

— 

4. 70 



6 

5. 05 

4.94 

5. 05 

4.88 

4.89 

6. 71 

— 

5. 10 



8 

5.03 

4. 92 

5. 27 

5, 13 

5.15 

7. 67 

— 



.. 

10 

4.88 

4. 74 

5. 42 

5. 32 

5. 34 

8.64 


5.58 

Ogive, 

3.00 

0 

2.69 

2.68 

2.68 

2. 60 

2.60 

2. 91 

2. 30 

2.70 

/d— 5 


2 

3.00 

3.00 

3. 22 

3. 03 

3. 01 

3.87 


3. 02 



4 

2.88 

2.88 

3.55 

3.28 

3.26 

4.85 


3. 35 



6 

2. 56 

2.56 

3.73 

3. 42 

3.39 

5.84 


3. 50 



8 

2. 30 

2.29 

3.83 

3. 49 

3. 46 

6. 82 





10 

2. 16 

2. 15 

3.88 

3. 53 

3.50 

7. 82 

— - 

3. 75 


4. 24 

0 

2. 74 

2.66 

2.66 

2.59 

2 59 

2. 84 

2.30 

2. 70 



2 

3.28 

3.21 

3.24 

3.08 

3.08 

3.74 


3. 15 



4 

3.48 

3.43 

3.62 

3. 42 

3. 43 

4.68 


3.50 



6 

3. 39 

3. 34 

3.87 

3.66 

3.67 

5. 65 


3.75 



8 

3. 12 

3.06 

4.03 

3.82 

3.82 

6.63 

— 


: 


10 

2.76 

2.70 

4. 13 

3. 92 

3.92 

7.61 


4.09 


5. 05 

0 



2.64 

2.58 

2. 57 

2. 79 

2.30 

2.70 



2 



3.23 

3.08 

3.09 

3. 66 


3.20 



4 



3.65 

3.48 

3. 49 

4.58 


3.73 



6 



3.94 

3.77 

3.78 

5.53 


3.85 



8 



4. 15 

3.98 

3. 99 

6.50 




__ 

10 




4.28 

4.12 

4. 14 

7. 47 


4.06 


6.28 

0 



2. 62 

2.55 

2.54 

2. 72 

2.30 

2. 62 



2 



3. 19 

3.05 

3.06 

3.54 


3.13 



4 



3.62 

3. 47 

3. 49 

4. 43 


3, 60 



6 



3.94 

3.80 

3.85 

5.36 


4.00 



8 



4.20 

4.06 

4. 12 

6.30 


— 


-- 

10 

---- 

— - 

4.34 

4. 26 

4. 33 

7.26 


4. 32 

Ogive, 

3.00 

0 



1.59 

1.56 

1.56 

1.70 

1.35 

1.65 

/□ — 3 


2 





2. 12 

2.03 

2.04 

2. 59 

— 

2.12 



4 




2. 43 

2. 34 

2. 35 

3.53 


2. 45 



6 



2.60 

2.52 

2. 54 

4.51 


2. 49 



8 



2.68 

2.62 

2.64 

5. 49 


— 


.. 

10 


— 

2. 73 

2. 67 

2. 69 

6. 47 


2.60 


4. 24 

0 



1. 56 

1.52 

1.52 

1.61 

1.35 

1.58 



2 



2. 11 

1.98 

2.01 

2. 43 



2.17 



4 


^ 

2.51 

2.34 

2.38 

3.34 



2. 44 



6 



2. 76 

2.58 

2.64 

4.28 


2.75 



8 



2. 93 

2. 74 

2.80 

5.24 




.. 

10 

— 

— 

3.03 

2. 84 

2. 91 

6. 22 


3.12 


5.05 

0 



1. 55 

1.51 

1.48 

1.58 

1.35 

1.47 



2 



2.26 

1.96 

1.97 

2.35 



2. 25 



4 



2.79 

2. 32 

2. 35 

3.23 



2.75 



6 




3. 19 

2. 59 

2.63 

4. 16 



3. 15 



8 




3. 49 

2. 77 

2.83 

5.11 







10 

— 

— 

3.73 

2.90 

296 

6,08 

.... 

3.60 


6.28 

0 



1.52 

1.50 

1.45 

1.53 

1.35 

1.53 



2 





2.06 

1.90 

1.88 

2.25 



2. 10 



4 





2. 45 

2.24 

2.23 

3.09 



2. 45 



6 




2. 71 

2.49 

2. 49 

3.99 

— 

2.66 



8 





2. 87 

2. 67 

2.67 

4.92 

— 

— 


-- 

10 

— 

.... 

2.98 

2. 79 

2.80 

i 5.87 

.... 

3.05 
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Positive directions of axes and angles (forces and moments) are shown by arrows 


Axis 

Force 

(parallel 

to axis) 
symbol 

I 

Moment about axis 

Angie 

Velocities 

Designation 

Sym- 

bol 

Designation 

Sym- 

bol 

■ 

Positive 

direction 

Designa- 

tion 

Sym- 

bol 

Linear 
(compo- 
nent along 
axis) 

Angular 

Longitudinal 

X 

X 

Rolling 

L 

Y >Z 

Roll 

ijp 

u 

V 

Lateral 

Y 

Y 

Pitching. ' 

M 

Z 

Pitch 

' e 

V 

q 

Normal 

Z 

Z 

Yawing ; 

N 

X ► Y 

Yaw 

4 

1 

1 w 

i 

r 


Absolute coefficients of moment Angle of set of control surface (relative to neutral 

n __ L p M ~ _ N position), 5. (Indicate surface by proper subscript.) 

Cn ~$s 

(rolling) (pitching) (yawing) 


4. PROPELLER SYMBOLS 


D 

P 

Diameter 
Geometric pitch 

P 

p 

Power, absolute coefficient 

P/D 

V’ 

Pitch ratio 
Inflow velocity 

c, 

5 ipv b 

Speed-power coefficient 

V, 

Slipstream velocity 

V 

Efficiency 

T 

Q 

» T 

Thrust, absolute coefficient a 

Torque, absolute coefficient 

n 

* 

Revolutions per second, rps 
Effective helix angle— tan -1 


5. NUMERICAL RELATIONS 

1 hp=76.04 kg-m/s=550 ft-lb/sec 1 lb=0.4536 kg 

1 metric horsepower — 0.9863 hp 1 kg=2.2046 lb 

1 mph=0.4470 mps 1 mi— 1,609.35 m — 5,280 ft 

1 mps=2.2369 mph 1 m = 3.2808 ft 








